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ABSTRACT 

In this paper we investigate statistical equilibrium of Ti in the atmospheres of late-type stars. 
The Ti I/Ti II level populations are computed with available experimental atomic data, ex- 
cept for photoionization and collision induced transition rates, for which we have to rely on 
theoretical approximations. For the Sun, the NLTE line formation with adjusted H I inelas- 
tic collision rates and MAFAGS-OS model atmosphere solve the long-standing discrepancy 
between Ti I and Ti II lines. The NLTE abundances determined from both ionization stages 
agree within 0.01 dex with each other and with the Ti abundance in C I meteorites. The Ti 
NLTE model does not perform similarly well for the metal-poor stars, overestimating NLTE 
effects in the atmospheres of dwarfs, but underestimating overionization for giants. Investigat- 
ing different sources of errors, we find that only [Ti/Fe] ratios based on Ti II and Fe II lines can 
be safely used in studies of Galactic chemical evolution. To avoid spurious abundance trends 
with metallicity and dwarf/giant discrepancies, it is strongly recommended to disregard Ti I 
lines in abundance analyses, as well as in determination of surface gravities. 

Key words: Radiative transfer - Line: profiles - Line: formation - Sun: abundances - Stars: 
abundances 



X 
S3 



1 INTRODUCTION 

Most of abundance analyses in metal-poor stars focus on few key 
elements, among others, Ti. This has several reasons. The main 
nucleosynthesis site of Ti has not yet been unambiguously iden- 
tified and all models of Galactic chemical evolution completely 
fail to describe observ ational trend of [Ti/Fe] with metallicity 
( lMcWilliametaI.ll 19951) , which increases with decreasing [Fe/H]. 
In spectra of late-type stars, Ti is represented by a large number of 
spectral lines in both neutral and ionized stages. Therefore, Ti can 
be used to verify basic stellar parameters (reff , log g) and micro- 
turbulence, which are determined, for example, by means of pho- 
tometric methods or from the excitation and ionization balance of 
Ffl 

Traditionally, Ti abundances in stellar atmospheres are de- 
termined assuming that distributions of atoms among excitation 
states and ionization stages derive from the formulas of Saha- 
Boltzmann for the local values of electron temperature and density 
in models of atmospheres, so-called local thermodynamic equilib- 
rium (LTE). Important evidence against this assumption comes, in 
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' In this context the notation of excitation and ionization equihbrium refers 
to the equality of abundances determined from all detected spectral lines 
of an element, which is represented by several ionization stages in a stellar 
atmosphere 



return, from LTE analyses of high-resolution spectra. It is known 
that low-excitation Ti I lines deliver systematically lower abun- 
dances compared to the high-excitation lines and lines of Ti II 
for di ff^erent stellar paramet ers. LTE studies of the Sun with ID 
static telackwell et al."l987*) and 3D hydrodynamical model atmo- 
spheres ( Asplund et al. 2009) demonstrate that these discrepancies 
are larger than 0. 1 dex, which is unacceptable given the precisely 
known fundamental parameters of the Sun. The same problem was 
reported f or other nearby sta rs with well-known parameters, such 
as Pollux jRulandetalJ 19801) . 

Analyses of metal-poor stars also point to pronounced NLTE 
effects on Ti I in their atmospheres. Systematically higher [Ti/Fe] 
abundance ratios in metal-poor dwarfs compared to giants with 
similar metallicites were reported by Bonifacio et al. (2009). Posi- 
tive offsets between abundances based on Ti II and Ti I lin es were 
reporte d for horizontal branch stars dClementini et al.l 1 995*) and for 
giants terownet alJI 19831: lGratton & SnedeiJll99ll :|johnson 200 J; 
'Tafelmever et al.ll2010l) . iLai et al.l (12008*) found a trend of [Ti/Fe] 
abunda nces ratios w i th ste ll ar T^g, which was supp orted by the 
data o f IPreston et alj ( l2006l) . ICavrel et al.l ( [2004h , and lCohen et"ai] 
( |2004|) . 

Theoretical studies of NLTE effects in Ti are very sparse. 
iHauschildt. Allard & Biroril h999l) investigated statistical equilib- 
rium of Ti in M-type dwarfs and giants. They found small effect 
of Ti I overionization on the number density of TiO molecules, 
the main opacity agent in the atmospheres of M stars. However, 
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they demonstrated that the NLTE effects on the Ti I fine formation 
are important and grow wit h increasing model Teff for dwarfs and 
decrea sing Teff for giants. Manso Sainz & Landi Degl'Innoc entil 
( l2002h and IShchukin a & Trujillo Bueno (2009) used NLTE pop- 
ulations of Ti to model polarization in the Ti I lines of multiplet 
42. Their goal was to study weak solar magnetic fields. None of 
the investigations focused on the NLTE effects on Ti abundances in 
stellar atmospheres. 

We present first results of a study aimed at the NLTE mod- 
elling of Ti line formation in the atmospheres of late-type stars. In 
Sect.|2l we give a detailed description of the methods developed for 
NLTE and spectrum synthesis calculations. Statistical equilibrium 
of Ti for a restricted range of stellar parameters is discussed in Sect. 
13.11 Ti abundances for the Sun and four metal-poor stars are given 
in Sect. [3] The results are summarized in Sect.|4] 



2 METHODS 

In order to carry out the abundance analysis, we use a unique pro- 
gram package developed at the University Observatory Munich 
(the group of T. Gehren). The package consists of three codes, 
which have been developed consistently and are adapted for spec- 
troscopic analysis of A-,F- and G-type stars. MAFAGS is a ID 
LTE model atmosphere code. Detail is a code for solving multi- 
level non-LTE radiative transfer problems with a given static ID 
model atmosphere. SIU is an interactive non-LTE spectrum syn- 
thesis code that makes use of MAFAGS model atmospheres and 
level populations from DETAIL. Below, we give some details for 
each of the programs. The general procedure of abundance de- 
ter mination is the same as in our previous work on Mn, Co, and 
Cr ([B ergemann & Gehren' 2008; Be rgemann. Pickering & GehrenI 
l20IQ : .Bergemann & CescuttL2010.) . 



2.1 Model atmospheres 

The static plane-parallel LTE model atmospheres were kindly gen- 
erated by F. Grupp with the MAFAGS code. In the most recent 
version, MAFAGS-OS, line blanketing is treated with an opacity 
sampling method, with line absorption sampled at ~ 86 000 wave- 
length points. All opacity sources relevant to modelling AO-KO 
stars including line absorption due to various diatomic molecules 
are taken into account. Full description of these models, includ- 
ing the reference set of solar element abundances, as well as 
selection crite ria of the wavelength grid and lines, is given in 
iGrupd ( 120043) . In particular, logepe = 7.50 was adopted, which 
is consistent with the NLTE solar abundance of Fe determined by 



Gehren et al. (2001). Convection is included in the formulation of 



Canuto & Mazzitellil ([T99ll) . which differs from the mixing-length 



theory of Bohm-Vitensd ( Il958l) in that it assumes a full spectrum 
of turbulent eddies. Compressible turbulence is accounted for us- 
ing the mixing length, A = (icm/Hp, where Hp is the pressure scale 
height and the free parameter a ^^ is set to 0.82. The latter value 
was derived bv lBemkopj i ll 9981) from stellar evolution calculations 
for the Sun and from the requirement that the MAFAGS models fit 
the observed Balmer line profiles. 

We also performed abundance calculations with the model 
atmospheres compute d with older version of the MAFAGS code 
dFuhrmann et aljl997 ): these models were used to determine spec- 
troscopic stellar parameters for the selected stars (Sect. 13. 3.21 ). 
This version is based on opacity distribution functions (ODE) from 
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Figure 1. Difference in temperature and gas pressure stratifications of 
the solar theoretical MAFAGS-OS, MAFAGS-ODF, and semi-emprirical 
MACKKL models. 



lKuruczl l fT993) . Convection was taken into account with the mixing- 
length theory ( lB6hm-Viteris3ll958l) and the mixing length was set 
to 0.5 pressure scale heig hts, the latter value was also calibr ated on 
the solar Balmer lines bv lFuhrmarm. Axer & Gehrei] l ll993h . A de- 
tailed comparison of the solar MAFAGS-ODF and MAFAGS-OS 
models is given in Grupp ( 2004b). 

MAGAFS models do not include chromospheres that is the 
case with all other theoretical ID models so far. Clearly, this 
is a crude approximation for the Sun. The existence of chro- 
mos pheres has also been demonstr ated for m etal-poor dwarfs 
(iSmith & Churchil l 1998) and giants jDupree. Hartmann & SmithI 
ll99Cl) . To discuss this deficiency of the models, we will also 
present some results for the semi-empiri cal solar model atmosphere 
with chromosphere jMaltbvetalJl 19861 . hereafter MACKKL). The 
MACCKL model was interpolated to 80 depth point s to increase 
numerical accuracy in line formation calculations dReetj [T9991 
PhD thesis). Differences between temperature and pressure strat- 
ifications of the solar MACCKL, MAFAGS-ODF, and MAFAGS- 
OS models as a function of continuum optical depth at 5000 A, 
logTsooo, are shown in Fig[T] 

We emphasize that the semi-empirical model atmosphere is 
used here only to demonstrate the shortcomings of the theoretical 
atmospheres, which are still the main tool to model any star other 
than the Sun. Thus, the solar abundance analysis, as well as the dif- 
ferential analysis of the metal-poor stars, are performed only with 
the MAFAGS models. 



2.2 Statistical equilibrium and line formation codes 

The NLTE level populations for Ti were computed with an updated 
version of the DETAIL code dButler & Giddingslll985b . In calcula- 
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tions of statistical equilibrium, each Ti line was treated with a Gaus- 
sian profile with 9 frequency points. The code solves a restricted 
NLTE problem, i.e. coupled statistical equilibrium and radiative 
transfer equations are solved for a fixed input model atmosphere. 
This may not be a bad approximation for Ti in the range of stellar 
parameters we are interested in. Ti I, which is, in fact, affected by 
NLTE, is not an important opacity source in the atmospheres of FG 
stars and does not contribute significantly to the free electron pool 
due to the low element abundance. Ti II is the dominant ionization 
stage and its number densitie s are well described in LTE (see Sect. 
[STT i. iHauschildt et alj (1 19991) computed full NLTE line-blanketed 
model atmospheres for M dwarfs and giants with solar metallic- 
ity including Ti as one of the NLTE species. They showed that the 
NLTE effects in Ti influence the atmospheres of cool stars only in 
the upper layers by changing the concentration of TiO+ molecules, 
and the effect on the number density of TiO molecules is negligible. 
Also, the effect of Ti on the atmospheric structure, even if present, 
is much smaller than that of Fe, which is an important source of the 
bound-free opacity in the UV and has, by far, the lar gest number 
of lin es all over the spectrum of a typical F-type star dGrupp et al.l 
llooi Fig. I). 

Still one has to keep in mind that, in general, the assumption 
of LTE in modeUing stellar atmosphe res is not realistic. As demon- 
strated bv lShort & Hauschildtl JioO^ . the integral effect of NLTE 
in Fe-group elements (Ti, Mn, Fe, Co, Ni) on the structure of the so- 
lar model atmosphere is significant. The differences in the tempera- 
ture structure of their NLTE PHOENIX model with respect to LTE 
models are as large as ±200 K at different depths. These differences 
are comparable with the temperature fluctuations i n the line form - 
ing regions in the 3D hydrodynamical solar model ( lAsDlunj2005h . 
Since NLTE effec ts in the presence of convectiv e inhomogeneities 
are amplified (Fe: IShchukina & Truiillo Buenoli2001) . one can ex- 
pect that self-consistent 3D NLTE hydrodynamical model atmo- 
spheres will have radically different structure compared to existing 
ID LTE, ID NLTE, or 3D LTE models. On the other side, such 
physically realistic models do not exist yet, and there is little chance 
they will appear in the near future. 

Emergent fl ux spectra w ere computed with an updated version 
of the code SIU l lReetzlll999h . using MAFAGS model atmospheres 
and departure coefficients of Ti I/Ti II levels from the DETAIL 
code. Ti lines were computed with full Voigt profiles taking into ac- 
count various broadening processes (see Sect. 13.2b . Ti abundances 
were determined by visually fitting the LTE and NLTE synthetic 
line profiles to the observed flux spectra including LTE modelling 
of blending features. This method is more reliable than ^-square 
fitting or equivalent width (EW) measurements for the solar-type 
stars, because a majority of Ti I and Ti II lines are blended and/or 
display an asymmetry in the core and red wing, which is a typical 
sign ature o f atmospheric temperature and velocity inhomogeneities 
(e.g. lAsplund et al. 2000) not taken into account by static plane- 
parallel model atmospheres. 



2.3 Ti model atom 

A model of the Ti atom was constructed with energy levels, wave- 
lengths of transitions, and transition probabilities from the Kurucz' 
databas^ which includes all laboratory data. The number of en- 
ergy levels is 216 for Ti I and 77 for Ti II, with uppermost excited 



levels located at 0.17 eV and 1.1 eV below the respective ioniza- 
tion limits, 6.82 eV and 13.58 eV. The model is closed by the Ti III 
ground state. The total number of radiatively-allowed transitions is 
4671 (3435 Ti I and 1236 Ti II). Fine structure was neglected in the 
statistical equilibrium calculations, except for the ground state of Ti 
I fl^F (configuration ls^2s-2p^3s^3p^3d^4s^). Excitation energy 
of each LS term was computed as a weighted mean of statistical 
weights and excitation energies of fine structure sub-levels. Also, 
transitions between fine structure levels were combined, with the 
total transition probability of each term being the weighted mean 
of loggf's of their fine structure components. A Grotrian diagram 
for Ti I is shown in Fig.[2l 

The electron collision cross-sections from states with allowed 
bound-bound and bound-free transitions were computed with the 
formulas of Ivan Regemorterl j 19621) and ISeatonI ( Il962h . respec- 
tively. Electron collision cross-sections from states connecte d only 
by f orbidden transiti ons were tre ated with the formula of I AllenI 
{19730. According to iMashonkinal j 19961) collision rates obtained 
with these formulas have only an order of magnitude accuracy. 

Threshold photoionization cross-sections for the levels y 
(3.93 eV) z^H" (3.96 eV), and v^F°(4.2 eV) were taken from 
lYang et al] ^200 91 and a hydrogenic approximation was assumed 
for their frequency variation. The cross-sections were measured 
by resonance ionization mass spectrometry and for the v ""F" levels 
they are con sistent with Har tree-Fock calculations with relativistic 
corrections /Saloman''l993). For the other levels, even quantum- 
mechanical calculations are not available, thus the photoionization 
cross-sections were derived from the hydrogenic approximation 
i MihalasI 1978t). T he experimental threshold cross-sections from 
lYang et alj l l2009h are lower than the hydrogenic cross-sections 
with elfective principal quantum numbers by one order of magni- 
tude. For example, the measured values for the v ^F° fine structure 
levels are 0.6 - 1.2 Mb, whereas the hydrogenic cross-section is 
roughly 18 Mb. 

In many aspects inclu ding the atomic s tructu re, our model 
atom is very similar to that of lHauschildt et all ( 1 19991) with the main 
differencq^ that inelastic collisions with H i atoms are included in 
our model. The corresponding boun d-bound and bo und-free rates 
were computed from the formulas o f lDrawirJ ( fT969l) in the version 
of Steenbock & Holweger ( 1984). 

Since the atomic data are of a low accuracy, we perform cal- 
culations of ionization equilibria for several Ti NLTE model atoms, 
which differ in the efficiency of e" and H I collisions (Sect. |T23] l. 
Thus, for the total e" collision rates we use the scaling factors 
5e = 0.01, 1, 10, and for the H I colhsion rates Sh = 0.05, 3. The 
final adopted values are 5h = 3 and 5e = 1, which give the smallest 
abundance scatter in the solar abundance analysis (see Sect. l3.2t . 



2.4 Line parameters 

A careful selection of Ti lines is essential, because ionization equi- 
librium of Ti in the solar photosphere is used to calibrate the effi- 
ciency of inelastic collisions with H I. 

Ti I and Ti II lines were s elected by the criti cal inspection of 
the KPNO atlas of solar fluxes ( Kurucz et al. "1984') and of the disk- 
center intensity spectrum (Brault & Testerman 1972). For Ti I, we 
rejected all lines that showed obvious blends in the profile, i.e. the 
contribution of the theoretically computed blend to the total Wa of 



^ http://kurucz.harvard.edu/atoms.html 



^ IHauschildt et al.l j 19991) also used different photoionization cross-section 
for the Ti I ground state. 
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a line was more than 5 %. For Ti II, this threshold was somewhat 
higher, because of a very small number of sufficiently strong Ti II 
lines with reliable transition probabilities. 

Some regions in the observed solar spectrum, such as shown 
in Fig.|3] suffer from excessive line blanketing and are not well re- 
produced by our model atmospheres and line lists. All Ti lines lo- 
cated in such regions were rejected from the solar analys is. For the 
other r egions, continuum is poorly defined. Following Ge hren et al.l 
( l200lh . we renormalized the continuum in these regions if the de- 
pression extended over few tenths of A. The re-adjustment of the 
continuum was never larger than ~ 0.5%, which is within the con- 
tinuum placement uncertainty. 

Tables[T]and|2]contain important information for each line in- 
cluding our evaluation of the line 'quality' in the solar spectrum, 
i.e. whether the line is blended or has an asymmetric profile. Term 
designations, wavelengths and excitation potentials of lower lev- 
els were taken fro m the NIST database, whenever available (oth- 
erwise from Martin . Fuhr & W iese 1988). The equivalent widths 
given were computed from the best-fitting synthetic profiles, ex- 
cluding contribution of blends. Equivalent widths are not used in 
the abundance calculations and serve only as a demonstration of 
line strength. For example, lines with Wa < 50 mA are not de- 
tectable in our spectra of metal-poor stars. 

2.5 Isotopic shift 

Ti is represented by five stable isotopes ''*Ti:'*^Ti:''^Ti:'"Ti:^''Ti and 
their relative abundances in the solar system are 8 : 7.43 : 73.8 : 
5.5 : 5.43 respectively. It has been recognized long ago that Ti lines 
in the solar spectru m are subj ect to an additional broadening due to 
the isotopic effect ( iAblll 19521) . In analogue to hyperfine splitting 
isotopic shift of energy levels leads to splitting of spectral lines 
into several components. For lighter multi-electron atomic systems, 
such as Ti, the isotopic shift of energy levels is primarily due to the 
differences in nuclear masses of isotopes giving rise to normal and 
specific mass shifts. As shown below, non-negligible errors in the 

"* http://www.mst.gov/pml/data/handbook/index.cfm 
^ Hyperfine splitting also affects atomic energy levels of '^^Ti and '^'Ti, 
but the fraction of these odd-Z isotopes is only 13% of the total element 
abundance. Thus, hyperfine splitting is ignored in this work. 



solar Ti abundance may arise if isotopic structure is neglected in 
the line formation calculations. 

Isotopic shifts in the UV and b l ue Ti II transitions were 
recently ca l culat ed by ISerengut et al.l ( l2008h and measured by 
iNouri et^ ( |2010 \ For some Ti II transitions, the total separation 
of components due t o the isotopic effec t is not larger than 6v ~ 900 
MHz or 6A~ 6 mA jNouri et al.ll201(]|) . In contrast, very large iso- 
topic shifts have been observed for transitions from the unfilled 3d 
shell (e.g. 3dMs^ 3<f 4p, or 3d 3d- 4p). For t he Ti II line 

at 448 8 A, which we use in our abundance analysis, iNouri et aP 
1 I2OIOI) derive an isotope shift of ~ 2400 MHz with the total sepa- 
ration of components of ~ 4500 MHz. This corresponds to a line 
width of ~ 30 mA that is comparable with thermal line broaden- 
ing in the solar atmosphere. There are also few measurements of 
the isotopic shift in Ti I, which include several transitions from our 
sample. 

In this study, isotopic structure was taken into account for two 
Ti I and three Ti II lines, which show large isotopic splittings, i.e. 
the separation of c omponents is ~ 15-20 mA. For the Ti I line s, we 
used the data from lGangrskv. Marinova & ZemlvanoH il99^ . The 
wavelengths and g/- values of individual components are given in 
Table [3] The total oscillator strengths for the transitions are given 
in Tables[T]and|2] 



3 RESULTS 

In the following sections, we present the results of statistical equi- 
librium calculations for Ti for a restricted range of stellar parame- 
ters and discuss whether the LTE approximation can be adopted for 
different combinations of T^g, logg, and [Fe/H]. We also compute 
NLTE abundances of Ti for the Sun and four metal-poor stars (Sect. 
13.3.21 ) and attempt to calibrate the poorly-known cross-se ctions to 
inelastic H I and electron collisions following the classical ' Prawinl 
( 1968) formalism to satisfy ionization equilibrium of Ti I/Ti II. De- 
viations from LTE and their effect on the abundances are discussed 
in terms of NLTE abundance corrections defined as: 
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Table 1. Parameters of the Ti I lines used for the spectrum synthesis. The abundances in columns 10 - 13 are given for ^/-values from Blackwell- Whitehead et 
al. (2006), if the latter are specified in column 7 (otherwise for the g/-values measured by the Oxford group, which are given in column 8, and exact references 
to the latter are given in column 9). 
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2 


b 


o^Fj 


z^f; 


0.00 


-1.07 


-1.06 


1 


79 


4.78 


4.82 


4.88 


4.91 


5192.98 


2 


b 


o^Fj 


z^F? 


0.02 


-0.96 


-0.95 


1 


85 


4.78 


4.82 


4.88 


4.91 


5210.39 


2 


b 


a^F^ 


z^F° 
' 4 


0.05 


-0.85 


-0.82 


1 


88 


4.84 


4.88 


4.94 


4.97 


5219.7 


2 


b 


a^Fj 


z^F" 


0.02 


-2.26 


-2.23 


1 


28 


4.88 


4.92 


4.97 


5.01 


5295.78 


31 


b 


O^Pj 




1.07 




-1.577 


3 


13 


4.82 


4.86 


4.92 


4.95 


5490.15 


45 


b 






1.46 




-0.877 


4 


24 


4.8 


4.82 


4.89 


4.91 


5662.15 


78 


b 


^ 4 


% 5 


2.32 




-0.053 


4 


24 


4.81 


4.83 


4.9 


4.93 


5866.46 


30 


b 


a^Pj 


y 3 


1.07 




-0.784 


2 


48 


4.81 


4.82 


4.9 


4.91 


5922.12 


27 


b,uc 


a\ 




1.05 




-1.410 


2 


21 


4.84 


4.87 


4.93 


4.96 


5965.83 


62 


b,uc 


a 




1.88 




-0.353 


4 


28 


4.8 


4.82 


4.89 


4.91 


5978.54 


62 


a 




z 


1.87 




-0.440 


4 


25 


4.79 


4.81 


4.87 


4.89 


6064.63 


25 


b,uc 


a^P 
u rjj 


z 


1.05 




-1.888 


3 


9 


4.85 


4.88 


4.96 


5.00 


6092.8 


61 


uc 




w^G" 


1.89 




-1.323 


4 


5 


4.82 


4.86 


4.92 


4.96 


6126.22 


25 


a 


a^P^ 


z^SJ 


1.07 




-1.369 


3 


23 


4.85 


4.87 


4.95 


4.97 


6258.71 


43 


a 


b^F^ 


.r 


1.46 




-0.299 


4 


53 


4.79 


4.79 


4.88 


4.89 


6261.1 


43 


b,uc 






1.43 




-0.423 


4 


49 


4.83 


4.84 


4.92 


4.93 


6303.76 


43 


b,rc 






1.44 




-1.510 


4 


9 


4.83 


4.86 


4.93 


4.97 


6312.24 


43 


b,rc 


b^F^ 


y'Gi 


1.46 




-1.496 


4 


9 


4.81 


4.84 


4.9 


4.94 


6554.24 


41 


b 


fo^Fj 




1.44 




-1.162 


4 


11 


4.83 


4.86 


4.91 


4.94 


6556.08 


41 


a 


b^F^ 


a-3f° 


1.46 




-1.018 


4 


11 


4.81 


4.83 


4.89 


4.91 


7357.74 


35 


b,rc 


fo^Fj 


y^^F^ 


1.44 




-1.066 


4 


23 


4.85 


4.87 


4.94 


4.96 


8426.51 


7 


b 


a'Fj 


z^^D° 


0.83 




-1.197 


2 


51 


4.89 


4.92 


4.96 


5 


8435.66 


7 


b 


a^F^ 


z'D° 


0.84 




-0.967 


2 


62 


4.88 


4.91 


4.96 


4.98 


8682.99 


24 


uc 


fl^Pj 


z^D^ 


1.05 


-1.69 


-1.82 


3 


15 


4.74 


4.77 


4.83 


4.86 



" Line identificati on: a - unblended; b - b lended; uc - uncertain continu um; s - saturated; rc - renorma lized continuum * Reference s to ^/-values: (1) 
iBlackwell etal] il982l a); (2) lBlackwell et alj h982l b); (3) lBlackwell et alj h983l) ; (4) lBlackwell et alj h986l) 
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Table 2. Parameters of Ti 11 lines used for the spectrum synthesis and individual abundances. The lines at 4394, 4395, 4443 A were not used in the solar 
abundance analysis. For some transitions iPickering et aL ( ,20Ql,) do not provide accuracies of g/'-values. 



A 


Mult. 


Rem." 


Lower 


Upper 




Wa 


log?/ 


Unc. 


Source'' MAFAGS-ODF 


MAFAGS-OS 








level 


level 


[eV] 


[mA] 






log ELTE 


log fNLTE 


log ELTE 


log ENLTE 


4394.05 


12 


b 






1.22 


77 


-1.78 


0.04 










4395.85 


3 


b 






1.24 


71 


-1.93 


0.04 


\ 








4443.8 


3 


b 






1.08 


138 


-0.72 


0.02 










4444.56 


6 


b 






1.12 


61 


-2.24 




1 4.97 


4.97 


4.99 


4.99 


4468.5 


6 


b 






1.13 


139 


-0.6 


0.02 


2 4.94 


4.92 


4.98 


4.96 


4488.32 


26 


b,uc 




x^F° 
^ ^3.5 


3.12 


50 


-0.51 


0.03 


1 4.86 


4.84 


4.91 


4.90 


4493.52 


2 


b 






1.08 


33 


-2.83 


0.05 


1 4.89 


4.89 


4.90 


4.90 


4563.77 


11 


b 






1.22 


129 


-0.69 


0.02 


1 4.98 


4.98 


5.01 


5.01 


4583.41 


8 


b 


a*P 




1.17 


31 


-2.92 




1 5.03 


5.03 


5.04 


5.03 


4636.32 


7 


b 


a*P 


,4po 


1.16 


20 


-3.02 




1 4.81 


4.8 


4.85 


4.84 


4657.2 


13 


b 


h%. 


,2po 


1.24 


52 


-2.32 


0.03 


2 4.98 


4.97 


4.99 


4.98 


4708.66 


10 


uc 




Z'F2.5 


1.24 


51 


-2.34 


0.05 


1 4.96 


4.95 


4.97 


4.96 


4798.53 


1 


a 




^'G2.5 


1.08 


44 


-2.68 


0.08 


1 4.97 


4.97 


5.00 


5.00 


5129.15 


20 


b 






1.89 


72 


-1.24 


0.03 


1 4.96 


4.93 


4.99 


4.96 


5336.81 


16 


b 




^'F3.5 


1.58 


71 


-1.63 


0.03 


2 5.05 


5.02 


5.03 


5.00 


5381.01 


16 


b 






1.57 


57 


-1.92 


0.05 


1 4.95 


4.94 


4.95 


4.94 


5418.77 


16 


b 






1.58 


51 


-2 


0.05 


1 4.90 


4.89 


4.90 


4.89 



' Line identification: a - unblended; b - blended; uc - uncertain continuum * References to g/-values: (l) IPickering etal.l i200ll) ; f2^ lBizzarri il993h 



3.1 NLTE effects on atomic number densities and line 
formation 

The qualitative analysis of level departure coe^cienjj^ helps to un- 
derstand NLTE effects in excitation and ionization balance of Ti. 
Departure coefficients of selected atomic energy levels of Ti I com- 
puted for MAFAGS-ODF model atmosphere^ with different pre- 
scriptions for elastic coUisional rates are shown in Fig. [S] Fig. |4] 
shows relative number densities N(Ti I)/N(Ti II) and N(Ti III)/N(Ti 
II) for the solar model atmosphere computed under LTE (dashed 
lines) and NLTE (solid lines). 

In general, most of the Ti I levels in the atmospheres of late- 
type stars are underpopulated relative to their LTE number den- 
sities. The driving mechanism is overionization by non-local UV 
radiation field with mean intensity 7,, larger than the local Planck 
function By(T^). However, due to the complexity of Ti atom, other 
NLTE mechanisms also play a role and population of each level 
at each atmospheric depth is determined by its collisional and ra- 
diative interaction with hundreds of other atomic states and, indi- 
rectly, even with the levels of the next ionization stage, Ti II. Thus, 
clear isolation of interaction processes leading to NLTE population 
of each atomic level is not possible and we have chosen only few 
Ti I and Ti II levels to illustrate the typical NLTE effects. Transi- 
tions between these levels are strong and few of them, /i4758 A 
and /i4759 A (a^H x^W), were selected for the stellar abun- 
dance analysis. Fig.|5]shows the Ti I ground state a^F, even-parity 
metastable levels a^H (isiow = 2.2 eV), odd-parity levels f^F° and 
x^H" (£|o„ = 4.8 eV and 4.85 eV, respectively). Also, few Ti II 



* The departure coefficients are defined as the ratio of NLTE to LTE num- 
ber densities of atoms in a certain excited l evel i, b, = (defini- 
tion according to lwiibenea & Zwaanll973) 

^ The choice of the ODF or OS model has no impact on the conclusions 
drawn in this section 




-5 -4 -3 -2-10 1 
log X (5000 A) 



Figure 4. Relative number densities N(Ti I)/N(Ti II) and N(Ti III)/N(Ti II) 
for the solar model atmosphere with parameters 5780/4.44/0 and [Ti/H] 
= 4.9, computed under LTE and NLTE. 

levels, which give rise to the Ti II transitions visible in our stellar 
spectra, are shown. 

Fig. |5^ illustrates departure coefficients computed for the so- 
lar mod e l atmo sphere assuming a scaling factor 5h = 0.05 for the 
iDrawinI (Il968h cross-sections for collisions with H I atoms and 
non-scaled cross-sections for transitions due to electron collisions, 
Se = 1 (Sect. |2.3"t . LTE distribution functions for the Ti I are invalid 
already at logTsooo ~ +0.3, which is where the atmosphere be- 
comes transparent to UV continuum radiation. The super-thermal 
radiation field, > ByiT^), increases photoionization rates from 
the Ti I levels with ionization edges at relevant wavelengths, that 
is, at /I ~ 300 - 400 nm. Since this process is not compensated by 
the recombination rates, which are fixed by the local kinetic tem- 
perature Te, the disbalance in number densities of well-populated 
Ti I levels with excitation energy 2 - 3 eV is transferred to the 
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-4 -3 -2-10 1 -4 -3 -2-10 1 -4 -3 -2 -1 

log T (5000 A) log X (5000 A) log x (5000 A) 





(h) 4800/2.2/-2.4 - 



-3 -2 -1 

log X (5000 A) 




-3 -2 -1 
log X (5000 A) 




-3 -2 -1 
log X (5000 A) 



Figure 5. Departure coefficients for selected Ti I and Ti II levels computed for the MAFAGS-ODF (a)-(d) and MACKKL (e)-(f) solar model atmospheres with 
dilferent prescriptions for inelastic collision rates, (g) to (i): Same as (b), but the model atmospheres are taken from the grid wiffi stellar parameters specified 
in each panel. Level designations in all panels are the same as in (a). See text for details. 



other levels via collisions and strong UV radiative transitions. This 
is an interesting case, when collisions inhibit deviations from LTE 
for some levels, but amplify them for the others by thermalizing the 
relative populations of the levels with different degree of underpop- 
ulation. 

At -1.5 < logTsooo < in the solar case, optical depth in the 
strong UV transitions of Ti I is still larger than unity and photoion- 
ization cross-sections of the levels in the hydrogenic approxima- 
tion are a monotonic function of frequency, cTph-ion ~ l/v^^, thus 
underpopulation of all levels is very homogeneous. The weak Ti I 
lines, which are formed at these optical depths, are only affected 
by the change in opacity, which scales with the departure coeffi- 
cient of their lower level. Since respective /7,-factors are less than 
unity, NLTE model lines are weaker compared to LTE and require 
larger Ti abundance to fit the observed spectrum. A few examples 
of NLTE synthetic lines of Ti I computed with the MAFAGS-ODF 
solar model are given in Fig. |6{a,b), where they are compared to 
the observed profiles. The shapes of line profiles computed with 
ODE or OS models are very similar, thus profiles computed with 
the latter are not shown in the Figure. The LTE profiles reproduce 
the NLTE profiles after adjustment of the Ti abundance and macro- 
turbulence. 

Excitation balance in Ti I significantly deviates from thermal 
at logTsooo ^ -1-5 for the MAFAGS solar model, where photon 
losses in numerous near-UV and blue transitions influence the rel- 
ative populations of the levels. The majority of such transitions are 



resonance (multiplet numbers 1 to 37). For example, the Ti I ground 
state, a ^¥ is radiatively connected with the t^V° level with the ex- 
citation energy 4.8 eV (Fig.[5]a). Although there is an ample pump- 
ing of the upper level by super-thermal /y at -1 < logTsooo < 0, 
photon escape in the lines of the multiplet 30 leads to a sudden un- 
derpopulation of /'F° at logTsooo 1-3. The latter process, com- 
bined with photon losses in other resonance transitions, causes the 
outward flattening of the departure coefficient for the Ti I ground 
state a^F, up to the depths where its fc, -factor drops again due to 
the dominant overionization. The same mechanism acts between 
the fl^H and x'H° levels (multiplet 114). The corresponding lines 
form at logTsooo ~ -1-2 (Fig.[5li). 

In the outermost layers, logTsooo ^ -3.5, the theoretical so- 
lar model atmosphere is not adequate: it lacks chromospheres. In 
the MACKKL semi-empirical solar model, the temperature mini- 
mum occurs at logTsooo ~ -3.4. Below this optical depth, the be- 
haviour of departure coefficients is identical to the MAFAGS model 
(Fig.|5^,f). However, the outward rise of temperature leads to the 
overpopulation of Ti I levels, contrary to the strong underpopula- 
tion predicted by the MAFAGS models. As a result, NLTE abun- 
dance corrections change in sign and magnitude, which is partic- 
ularly important for lines formed above T^-^\^, such as the strong 
resonance Ti I lines at 3998.6 A and 3989.76 A. LTE abundances 
determined with the MACKKL and MAFAGS models also differ 
because the former is slightly warmer than the latter in the layers 
below logTsooo = -3 (Fig.[T). 
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4443.6 4443.7 4443.8 4443.9 4444.0 4395.7 4395.8 4395.9 4396.0 

X(k) X(k) 



Figure 6. Synthetic (black trace) and observed (symbols) profiles of selected Ti I and Ti II lines in the solar spectrum, (d): for demonstration purposes, the 
theoretical profile was synthesized excluding the contribution of blends in the wings that makes the fit worse, (e): the line was not used for the solar abundance 
calculations due to an unknown blend. See text for details. 



To illustrate the difference between the theoretical and semi- 
empirical models with respect to abundance determinations, we 
consider two Ti 1 lines. The line at 4981.73 A (multiplet 44) is 
the strongest in our list with W,i ~ 120 mA. The NLTE line profiles 
computed with the MAFAGS-ODF and MACKKL solar model at- 
mospheres are compared with the observed solar flux profiles in 
the Fig.|6};. The NLTE abundance corrections have different sings: 
Anlte(MAFAGS - ODF) = +0.02 dex and Anlte(MACKKL) = 
-0.05 dex. The NLTE Ti abundances deduced from such profile fits 
are 4.86 dex and 4.98 dex for the MAFAGS-ODF and MACCKL 
models, respectively. For the weak line at 4997 A with Wa ~ 34 
mA, MACKKL model also overestimates the abundance by 0. 1 dex 
compared to the MAFAGS-ODF model, but the NLTE corrections 
for both model atmospheres are equal, Anlte = +0.04 dex. 

The influence of uncertainties in atomic data on the excita- 
tion and ionization balance of Ti was tested by changing the cross- 
sections for collision-induced transitions by hydrogen atoms and 
free electrons. Scaling down the cross-sections of transitions due to 



e" collisions by two orders of magnitude (Fig. [5};), S^. = 0.01, leads 
to decoupling of all levels at logrsooo > -1.5, whereas increas- 
ing the cross-sections by a factor of ten (Fig.|5}l) produces relative 
thermalization of the levels with excitation energy below 5 eV in 
the same range of optical depths logTsooo- 

Departure coefficients are very sensitive to the variation of the 
scaling factor to inelastic H I collisions 5h. With 5h = 3, the Ti 
I levels nearly thermalize at -2 < logTsooo ^ (Fig. [Sj)). It is 
only in the uppermost layers with low densities, logTsooo ^ -3, 
that the influence of collisions with H I is small. We note that 5h = 
3 produces the smallest abundance scatter between Ti I lines and 
simultaneously satisfies ionization balance of Ti in the solar case 
(Sect. |T2t . Thus, the latter value defines our reference model atom. 
Yet, it has to be kept in mind that calibration of 5h on the observed 
stellar spectrum may hide or compensate other deficiencies in the 
modelling, such as those related to the use of the ID static model 
atmospheres with a certain prescription for convection, i.e. mixing- 
length theory in our case, and the fixed microturbulence parameter. 
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Table 3. Isotopic components of selected Ti 11 and Ti I transitions for the 
solar mixture of the isotopes. See text. 



Transition Wavelength logg/ Isotope 

A 



Till 



4488.346 


-1.607 


46 


4488.337 


-1.638 


47 


4488.330 


-0.642 


48 


4488.322 


-1.77 


49 


4488.315 


-1.78 


50 



4468.498 


-1.697 


46 


4468.504 


-1.729 


47 


4468.510 


-0.732 


48 


4468.515 


-1.860 


49 


4468.521 


-1.868 


50 



Til 



4444.548 


-3.337 


46 


4444.554 


-3.369 


47 


4444.560 


-2.372 


48 


4444.565 


-3.500 


49 


4444.571 


-3.507 


49 


5866.440 


-2.108 


46 


5866.449 


-2.100 


47 


5866.459 


-0.972 


48 


5866.469 


-1.977 


49 


5866.480 


-1.937 


49 



5922.100 


-2.734 


46 


5922.109 


-2.726 


47 


5922.119 


-1.598 


48 


5922.129 


-2.603 


49 


5922.140 


-2.563 


49 



diative interaction between levels, i.e. non-balanced excitations at 
the depths where photons escape in the wings of strong Ti I lines 
and spontaneous de-excitations at the core formation depths. In 
contrast, very homogeneous distribution of /7,-factors with optical 
depth is characteristic of warm high-gravity models (Fig.|5j). This 
is an expected result, because collisions provide close coupling of 
the levels and even the Ti I ground state with the largest occupation 
number is subject to overionization caused by strong non-local UV 
radiation field near its ionization threshold. 

Deviations from LTE in Ti II are small for the combinations of 
stellar parameters investigated in this work. Departure coefficients 
of low and intermediate-excitation Ti II levels deviate from unity 
in the upper atmospheric layers, at log T5000 ^ -2. Most of the Ti II 
lines form below this depth, where fc, -factors are slightly larger than 
unity, and the corresponding effect on abundances is never larger 
than -0.05 dex. 

3.2 The Sun 

3.2.1 Spectrum synthesis 

The solar abundance of Ti was determined by visually fitting syn- 
thetic spectral lines computed with LTE or NLT E number densi- 
ties to the Solar Flux Atlas of iKurucz et alj ( 1 1984.1 . We used both 
MAFAGS-ODF and MAFAGS-OS model atmospheres. The stan- 
dard line broadening mechanisms were taken into account: so- 
lar rotation with Koi.o = 1-8 km s"', microturbulence velocity 
= 0.9 km s"', and a radial-tangential macroturbulence velocity 
^Rx = 2.5 ... 4 km s"' . Line broadening caused by elastic collisions 
with H I atoms was compu ted using the velocity pa rameters and 
temperature exponents from lAnstee & O'Maral ( Il995l) . In Tables □ 
and|2l this parameter is given in terms of the van der Waals damp- 
ing constant log Ce for each transition. 



Thus, our final value of 5h may not be appropriate for other NLTE 
studies of Ti, which are based on different model atmospheres and 
methods (see e.g. discussion in Asplund 2005). 

We performed test NLTE calculations for several models with 
the following parameters: (a) Teff = 4800, iogg = 2.2, [Fe/H] 
= +0.6; (b) Teff = 4800, logg = 2.2, [Fe/H] = -2.4; (c) T^g = 
6200, logg = 4.6, [Fe/H] = -2.4. These parameters are represen- 
tative of stars commonly used in Galactic chemical evolution stud- 
ies, i.e. metal-rich A^-type giants in the bulge, metal-poor giants 
and dwarfs in the halo. Departure coefficients for the three models 
computed with Su = 3 and = 3 are shown in Fig. |5] A rather 
counter-intuitive result is that the Ti I and Ti II level populations 
show departures from LTE even for the cool giant with super-solar 
metallicity (Fig.|5^). Even though UV radiation field is weakened 
due to low photospheric temperatures and excessive line blanket- 
ing, collisions at log g ~ 2 are too weak to ensure LTE conditions. 
Since the stellar flux maximum is shifted to longer wavelengths, 
ionization from the Ti I levels with higher excitation potential, e.g. 
a ^H, becomes more important compared to the solar case. Qualita- 
tively, our results for cool stars of solar metallicity are very similar 
to|Hauschildt et al. ( 1999, their Fig. 4). 

NLTE effects are particularly large for the metal-poor models 
with low gravity (Fig.|5ji). NLTE number densities of the majority 

of Ti I levels are reduced by a factor of 2 already at log T5000 0-5. 

The undulations of the ii,-curves with log T5000 are due to strong ra- 



3.2.2 Oscillator strengths 

Oscillator strengths for the Ti II transitions were adopted from 
IPickering. Thome & PerezI feOOll) . who combined their branching 
ratios obtained by the methods of Fourier transform spectrometry 
with the radiative lifetimes measured by means o f time-resolved 
laser-induced fluorescence bv lBizzarri et alj jl993h . The accuracy 
of the latter i s abo ut 5%. We also used the logg/ values of 
iBizzarri et al. I l ll993h . The oscillator strengt hs from both sources 
agree well with each other for log g/ > - 1 .5 jPickering et alj200ll 
their Fig. 3), but significant discrepancies, up to 0.3 dex, are present 
for the weaker transitions. The uncertainties of g/-values from both 
references are typically less than 0.05 dex, but for some transitions 
experimental error was as large as 0.2 dex. 

For the Ti I transitions, we relie d on the oscillator strengths 
of iBlackwell- Whitehead erZH2006l), wh o appl ied the same ex- 
perimental techniques as Ipickering et al.l j200lh . As a test case, 
we also used the g/-values of Blackwell et al. (1982a, 1982b, 
1983, 1986), which we correc ted by +0.056 as recommended by 
iGrevesse. Blackwell & Petfor j ([1989), who renormalized the rela- 
tive oscillator strengt hs measured with the Oxf ord furnace to more 
accurate lifetimes o f iRudolph & Helbid imj) . The estimated er- 
rors of the absolute /-values are less than 5%. 

All parameters of the selected lines, including diff^erent sets of 
oscillator strengths, are given in Tables [T] and |2] Given the above- 
mentioned discrepancies between oscillator strengths from differ- 
ent experimental sources, we decided to compute the average Ti 
abundance for each set of g/-values. 
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Table 4. Solar abundance of Ti computed with different NLTE model atoms and model atmospheres. The brackets specify a combination (Sn.Se) used for 
each model atom. See texts. 



Ion 


N 


gf" 


MAFAGS-ODF 
LTE (3, 1) 


LTE 


(3,0.01) 


MAFAGS-OS 

(3,1) 


(3, 10) 


(0.05,0.01) 


(0.05, 1) 


(0.05, 10) 


Til 


12 


1 


4.80 ± 0.05 


4.84 ± 0.05 


4.90 ± 0.05 


5.01 ±0.08 


4.94 ± 0.05 


4.95 ± 0.05 


5.11 ±0.06 


5.02 ± 0.05 


5.01 ± 0.05 




52 


2 


4.81 ±0.03 


4.84 ±0.03 


4.91 ± 0.04 


4.95 ± 0.06 


4.93 ± 0.04 


4.93 ± 0.05 


5.09 ± 0.09 


5.00 ± 0.05 


5.00 ± 0.04 


Till 


10 


3 


4.94 ± 0.06 


4.93 ± 0.07 


4.96 ± 0.06 




4.95 ± 0.06 












11 


4 


4.98 ± 0.04 


4.97 ± 0.04 


4.99 ± 0.04 




4.98 ± 0.04 











° References: (l) lBlackwell-Whitehead et alj <20()6l) : ( 2) Blackwell et al. (198 2a. 1982b. 1983, 1986). s caled by +0.056 dex following l&evesse et al.l 

i l 19891) : (3) |Pickering etalj i200ll) : (4) lBizzarri efZI il993h 



3.2.3 Abundance ofTi in the solar photosphere 

The mean Ti abundances obtained under LTE and NLTE with 
different model atmospheres, model atoms, and sets of oscillator 
strengths are given in Table |4] The individual abundances for each 
spectral line are listed in Tables [T] and |2] NLTE abundances given 
in these tables were derived assuming efficient inelastic H - Ti col- 
lisions in statistical equilibrium calculations, 5h = 3. The choice of 
the scaling factor will be described below. 

In Table ID it can be seen that the LTE abundances deter- 
mined from the Ti I lines are systematically lower than the NLTE 
abundances, with the difference of 0.05 - 0. 1 dex, which depends 
on the treatment of inelastic collisions in statistical equilibrium 
calculations. Also, the MAFAGS-ODF model atmosphere under- 
estimates the solar Ti abundance by ~ 0.1 dex compared to the 
MAFAGS-OS model. This is due to differences in the temperature 
and pressure structure of the models (Fig .[T|l. The NLTE values de- 
termin ed from the Ti I lines using the Blackwell- Whitehead et al.l 
(l2006h and Blackwell et al. (1982a, 1982b, 1983, 1986) gf values 
agree well with each other. Mean LTE abundances deduced with 
eithe r of the gf sets are also consistent. ^lackwell-Whitehead et al. 
( I2OO6 ) noted tha t the ^/-values for sev eral near-IR transitions of 
Ti I measured bv lBlackwell et al] ( Il983h are inaccurate due to non- 
equilibrium effects in the Oxford furnace, and should be modified 
by ~ 50%. Our results do not allow to confirm or refute this state- 
ment, since in both cases the individual abundances for the near-IR 
transition a^Pj - z^D^ (8683 A) differ by more than Icr from the 
mean value computed for a corresponding logg/ set (Table[TJ. 

NLTE abundances for each Ti I line derived with the 
MAFAGS-OS model atmosphere and different atomic models (de- 
scribed in Sect. l2.3b are shown in Fig.[7]as a function of line equiva- 
lent width, oscillator strength, and excitation potential of the lower 
level of the transition. The combination of the scaling factors 5h 
and 5c used for each model atom is indicated in the first panel 
of each row. The case of (0.05, 10) is omitted because the results 
are almost identical to the case (0.05, 1) (Fig.|7^). LTE results are 
also presented. Filled symbols denote the values computed with the 
logg/ set from Blackwell et al. (1982a, 1982b, 1983, 1986), and 
open symbol s correspond to the individua l abund ances derived with 
the data from lBlackwell- Whitehead et al.1 j2006h . There are 12 lines 
in common between the two experimental studies. Dashed lines 
show the average NLTE abundance determined for each atomic 
model with the g/- values of Blackwell et al., Q. We do not present 
analogous plots for the Ti Il-based abundances because NLTE ef- 
fects on Ti II in the solar case are negligibly small, which is evident 
from Table |4] 

Inspection of NLTE logepi values in Fig. |7] shows that NLTE 
abundances are very sensitive to the atomic data used in statistical 



equilibrium calculations. For some combinations of 5h and 5e, e.g. 
(3,0.01)(Fig.|7}i) and (0.05,0.01) (Fig.|7}l) the scatter between the 
Ti I lines is very large. The individual abundances differ by a fac- 
tor of two, although there are no physically significant trends with 
the equivalent width. In addition, abundance seems to be correlated 
with the oscillator strength and excitation potential of the lower 
level, if only the g/-values of Blackwell et al. (1982a, 1982b, 1983, 
1986) are considered. The g/-log e and Et^^-log e slopes are larger 
for the cases (a) and (d) that clearly gives rise to the large spread of 
abundance with ly.i. On the one side, this is an irrefutable argument 
against small scaling factors to e" collision rates, e.g. 5e = 0.01. 
One could also argue that small 5h values are inappropriate, be- 
cause they lead to uncomfortably large solar Ti abundances (Ta- 
ble|4), > 5.00, which are not consistent with the meteoritic valu^ 
However, the latter depends on the solar abundance of the refer- 
ence element Si, which ma y be subject to usual problems of spec- 
troscopic methods as well l lAsDlunjr200(]| : IShi et alj|2008l) . A sec- 
ond possibility is that the Ti model atom is incomplete in terms of 
number of levels and transitions. Mashonkina et al. (2010) have re- 
cently demonstrated that it is important to include highly-excited 
predicted levels and transitions between them in statistical equilib- 
rium calculations for Fe. Although these data also exist for Ti, the 
number of predicted levels and transitions exceeds 14000 and few 
millions, respectively. Even more levels and radiative transitions 
are predicted for Fe I: ~ 37 500, respectively, ~ 6 millions. Such 
atomic models are not tractable even with ID NLTE codes and re- 
quire efficient algorithms to reduce them by combining the levels 
and lines. We are currently working on this. 

It has to be kept in mind that photoionization cross-sections 
are also not appropriate for Ti. However, in the absence of quantum- 
mechanical calculations, it is impossible to predict whether the hy- 
drogenic approximation over- or underestimates the cross-sections. 
We conclude this in analogy to the similar ato m Cr I , for which 
the calculated photoionization cross-sections of Nahaj ( 12009,) dif- 
fer from the hydrogenic approximation by orders of magnitude, 
being significantly larger for some levels and smaller for others. 
In principle, the same is true for the electron and H I collision 
rates. Comparison of R-matrix calculations for e" collisions for 
Ca II jMelendez. Baiitista & Badnelll l2007h with the formula of 
Ivan Regemortei ( 19621) also reveals order of magnitude differences. 
What concerns inelastic H I collision, compared to the Drawin's 
formulas, ab initio quantum-mechanical calculations predict signif- 
icantly lower collision rates for certain transitions of s imple alkali 
atoms teelvaev & Barklemll2003l : iBarklem et al.1120101) . and they 



log£ = 4.91 ± 0.03 dex jAsplund et alj2009l) 
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Figure 7. NLTE Ti abundances computed with the MAFAGS-OS solar model and dilferent model atoms. Filled symbols denote the values computed 
with the logg/ set from Blackw ell et al. (1982a, 1982b, 1983, 1986), and open symbols correspond to the abundances derived with the data from 
iBlackwell-Whitehead etal] j2006l) . 
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show that, in addition to excitation, other effects hke ion-pair for- 
mation become important. 

There are also other explanations for the apparent trends of in- 
dividual abundances with the parameters of Ti I lines and/or large 
scatter, which are not related to the deficiencies of NLTE mod- 
elling. First, one can question the reliability of statistics based on 
small number of lines, i.e. there are just few lines in our list with 
small log^/'s. Also, the range of low-level excitation potentials 
is limited to - 2.5 eV. Second, the uncertainties of g/-values 
from Blackwell et al. (1982a, 1982b, 1983, 1986) might be under- 
estimated. It is interesting that the same regularity, i.e. variation 
of abundance with m ultiplet, was also foimd in the LTE analysis 
of solar Ti I lines by iBlackwell et al.l ( 1 19871 . their Fig. 4), and it 
is present for the individual Ti abundances based on the LTE as- 
sumption (Fig. 13) and/or the MAFAGS-ODF model atmosphere. 
In addition to uncertainties in the gf values measured with the Ox- 
ford furnace, Blackwell et al. ( 1987) also suggested that some Ti 
I lines are affected by photospheric temperature and velocity fluc- 
tuations, which are very pronounce d in 3D hydrodynamical solar 
model atmospheres ( lAsDlundll2003) . Although this is most likely 
the case for low-excitation temperature-sensitive lines, we do not 
see any trend of abundance with line strength (Fig.|7^). The lat- 
ter probably indicates that a single depth-independent microturbu- 
lence is not a bad approximation. Otherwise, the medium-strong 
lines {W,i ~ 60 . . . 100) being very sensitive to would show sys- 
tematically higher (or lower) abundances. 

The abundances derived from the NLTE and LTE profile fit- 
ting of the Ti II lines are given in Table |4] Although the standard 
deviations of the mean abundances, especially for the ^/-values of 
IPickering etd] ( 1200 ll) . are quite large, there is no trend of indi- 
vidual line abundances with E\„^, and logg/ (not shown). As 
expected from the results of Sect. l3.1l LTE abundances are slightly 
larger than that computed under NLTE. The effect of the isotopic 
shift on the abundance determinations is not negligible. For exam- 
ple, a single component used to fit the observed Ti II line at 4488 
A, overestimates the abundance by 0.06 dex compared to the re- 
sult obtained with 5 isotopic components (Fig.|6}l). To demonstrate 
that the profile composed of several components is also slightly 
asymmetric, the Ti II line was plotted excluding the contribution of 
blends in the wings. 

Some Ti II lines require very large abundances, log£ 5.05 
dex, to fit their profiles in the solar flux spectrum. Few of them 
also demonstrate broad asymmetric cores in the solar disk-center 
intensity spectrum, which can not be reproduced even taking iso- 
topic structure into account. O ne of these lines, 4443 A, with an 
isotopic shift of ~ 380 MHz ( iNouri et al.ll20IOl) is shown in Fig. 
[6^. Very specific profile shape s are probably not due to convective 
motions jAsplund et al]|2000l) . Thus, we neglect such lines in the 
solar analysis assuming that they contain unresolved blends in the 
inner core^. The synthetic profiles of other discrepant lines sys- 
tematically underestimate depths of both wings, e.g. 4395 A (Fig. 



' The presence of blends can be revealed either by searching for asym- 
metries in the observed line proiiles or by inspecting theoretical line lists. 
One problem is that often weak blends overlap with the core of a stronger 
line under investigation that leaves the profile of the latter fully symmetric. 
Another problem is that theoretical line lists are compiled from transitions 
observed in laboratory conditions, which do not correspond to astrophysi- 
cal conditions. A comparison of observed and synthetic spectra proves that 
many observed lines still lack identification, which may well contribute to 
the well-known "missing UV opacity" problem in the UV. It is natural to 
expect missing opacity in other spectral regions as well. 



[6f). Thus, the abundance estimate inferred from such profile fits 
is rather subjective and depends, among other parameters, on the 
adopted value for the macroturbulence, which otherwise would not 
affect determination of the abundance being an external line broad- 
ening mechanism. We refrain from using these lines in the solar 
analysis, also because they are very sensitive to 

In TableEl it can be seen that MAFAGS-OS model and NLTE 
line formation with 5h = 3 achieve good ionization equilibrium 
of Ti I/Ti II for the Sun. For the given 5h, the mean abundances 
determined with = 0.01, 1, 10 are also consistent with the Ti 
Il-based abundances within their respective uncertainties, however 
5e = 1 provides the smallest abundance spread for the solar case. 
With respect to the Ti I lines, almost equal mean abundances were 
obtained for the g/- values from Blackwell- Whitehead et alj ( |2006|) 
and Blackwell et al. (1982a, 1982b, 1983, 1986), the latter data 
scaled by -f-0.056 dex. However, our results for the Ti II lines 
favour the g/-values of Pickering et al. (2001), which lead to a bet- 
ter agreement between the two ionization states. The MAFAGS-OS 
NLTE abundances, whether based on Ti I or Ti II lines, also agree 
well with the Ti abundan ce in C I meteorites, log£ = 4.91 ± 0.03 
dex jAsplund et alj2009h . 



3.3 Ti ionization equilibrium for the metal-poor stars 

3.3. 1 Observations and stellar parameters 

To test ionization equilibrium of Ti in the atmospheres with re- 
d uced metal content, we have chosen four stars from the sample 
of'Bergemann & Gehrenl ( l2008h . The reader is referred to this paper 
for a detailed description of observations. The spectra of HD 84937, 
HD 14028 3, and HP 122 563 were taken from the UVESPOP sur- 
vey CBagnulo et al.ll2003h . The spectra of HD 102200 were taken 
by T. Gehren and colleagues on the ESO UVES echelle spectro- 
graph at the VLT UT2 in Paranal, Chile. The UVES spectra have a 
slit-determined resolution of /i/A/1 ~ 50 000 and a signal-to-noise 
ratio better than SjN ~ 300 near 5000 A. 

Stellar parameters were adopted from iGehren et alj j2004l 

20061) for HD 84937, HD 140283, and HD 102200, and from 



Mashonkina et al. (2008) for HD 122563. These studies used the 
same observed spectra, model atmospheres (MAFAGS-ODF), and 
the codes as adopted in this work that secures consistency of 
our abundance analysis. The effective temperatures Tea were de- 
termined from B aimer line profile fits, and gravities logg were 
based on Hipparcos parallaxes from the I eSAI j 19971 ) catalogue. 
The spec troscopic temperatures are consistent with Teff's deter- 
mined by Casagrande et al.l ( 120101) using the Infrared Flux method. 
ICasagrande et al.1 ( l2010 l) derived T^s = 6155 K for HD 102200, 
Teff = 5777 K for HD 140283, and T^ff = 6408 K for HD 
84937. The iron abundance and microturbulence velocities were 
derived from LTE fitting of Fe II lines, which do not suffer from 
NLTE effects for the ran ge of stellar parameters investigated here 
dMashonkina et alj|2010l) . 

MAFAGS-OS model atmospheres are available only for HD 
84937 and HD 122563 (L. Mashonkina, private communication). 
The gravities adopted for the MAFAGS-OS models are on average 
~ 0.1 dex larger than that used to compute the MAFAGS-ODF 
model atmospheres. The reason is that new reduction of Hipparcos 
data ( van Leeuwe n 20 7) rev ealed some errors in the earlier version 
of the catalogue ( IESA||1997|) . However, for our stars the parallax 
differences between both catalogue versions are within the parallax 
error itself, which propagates into 0. 1 dex error in log g. 

Stellar parameters are given in Table |5] The estimated errors 
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Table 5. Stellar parameters for the selected sample. See text. 
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MAFAGS-OS 




[ff/Fe] 






[Fe/H] 


Teff \ogg [Fe/H] 


km/s 
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6350 4.09 -2.15 
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0.4 
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Figure 8. NLTE (red symbols) and LTE (blue symbols) abundances of 
Ti in the metal-poor star s. For comparison, the LTE abundance ratios of 
lOratton & SnedeJ h99lh are shown with open and filled black symbols, 
respectively. 



are 100 K for T^ff, 0.1 dex for logg, 0.1 dex for [Fe/H], and 0.2 
km/s for f,. 



3.3.2 Ti abundances in the metal-poor stars 

The analysis of the metal-poor stars is differential with respect to 
the Sun. An individual abundance derived for each line in a stel- 
lar spectrum is related to the solar abundance from that line com- 
puted with the same parameters and assumptions in the spectrum 
synthesis, including model atmosphere type (OS or ODF). The Ti 
abundance in a star relative to the Sun is given by: 

[Ti/H] = log(g/£). - log(g/e)o 

We also used several Ti II hnes (4394.05, 4395.85, 4443.8 A), 
which have been excluded from the solar abundance analysis be- 
cause of profile distortions, which are likely caused by blends 
(Fig. |6^) and sensitivity of profile wings to microturbulence (Fig. 
|6f). However, in spectra of very metal-poor stars the lines are 
much weaker and appear to be symmetric showing no evidence for 
blends. Here, we discuss the stellar Ti abundances normalized to 
the Fe abundance of a star, [Ti/Fe]; this parameter is relevant to 
Galactic chemical evolution studies. 

The Ti abundances determined for four metal-poor stars are 
given in Table [6] Fig. [8] displays the difference between [Ti/Fe] 
ratios computed from the lines of two ionization stages, [Tin/Tii] 
= [Tin/Fe] - [Tij/Fe], as a function of stellar metallicity. NLTE 
and LTE ratios for the metal-poor stars are shown with red 
and blue symbols, respectively. The figure also includes LTE Ti 



abundances in metal-po or giants and dwarfs from the study of 
iGratton & SnedeiJ ( Il99 1[) . They define stars with log g < 3 giants. 

Our LTE abundances based on Ti I and Ti II lines agree within 
their respective uncertainties for three stars, except for the very 
metal-poor giant HD 122563. In contrast, ionization equilibrium 
is not satisfied under NLTE with 5h = 3 and 5e = 1. Except for 
the metal-poor giant, the Ti I lines give systematically higher abun- 
dances than Ti II lines. Assuming any other combination of colli- 
sion scaling factors makes the discrepancy for the three stars worse 
because NLTE effects on Ti I increase more for the metal-poor stars 
than for the Sun. Variation of Sc has almost no effect on the Ti abun- 
dances, whereas decreasing the elficiency of hydrogen collisions by 
few orders of magnitude, 5h = 0.05, leads to extreme overioniza- 
tion from the lowest Ti I states. This is expected because concen- 
tration of free electrons in metal-poor atmospheres is very low and 
thermalization of levels can only occur due to H I - atom collisions. 
As a result, the NLTE abundances determined with 5h = 0.05 are 
on average 0.2 dex larger than those computed with 5h = 3, and 
the difference between the latter and the LTE-based abundances is 
of the same order. 

For HD 84937, our Ti Il-based LTE abundance is 0.15 dex 
lower than the LTE abundance determined from the Ti I lines. This 
discrepancy has no relation to NLTE effects, because there is no 
physical mechanism to produce overpopulation of low-lying Ti I 
levels at expense of Ti II in the atmospheres, where Ti I is the minor 
ionization stage. For HD 122563, we encounter a different problem: 
the NLTE abundance determined from the Ti I lines is 0.2 dex lower 
than that from the Ti II, although the discrepancy is twice as small 
compared to LTE. 

Thus, it appears that our reference NLTE model overesti- 
mates NLTE elfects on Ti in high-gravity atmospheres, corre- 
sponding to dwarfs, however it underestimates overionization in 
the atmospheres of giants. For metal-poor dwarfs and subgiants, 
nearly thermalized occupation numbers for Ti I levels are neces- 
sa ry to satisfy ionization b alance. As seen in Fig. [8] the results 
of iGratton & SnedenI l [l99lh support this conclusion. In their LTE 
analysis of metal-poor stars, the mean difference between ionized 
and neutral species is +0.07 ± 0.03 dex for dwarfs and +0.2 ± 0.03 
dex for giants. As noted above (Sect.[T), similar offsets fo r metal- 
poor g i ants were also r eport ed in the LTE studies b y .Brown et al.l 
i ll983l) , |johnsoril ( l2002h . and lTafelmever et al.l |201(]|) . 

3.3.3 Discussion 

The analysis of four metal-poor stars showed that the LTE approxi- 
mation in Ti line formation calculations is not adequate for evolved 
stars. However our NLTE models with available atomic data for Ti 
do not perform better if both evolved and unevolved stars are con- 
sidered. In fact, there are few reasons why the problem can still be 
attributed to deficiencies of the NLTE model. 

All lines detected in the stellar spectra originate from the 
metastable levels a^F^^j with the excitation energy ~ 0.8 eV. 
These levels, among many other low-lying states, dominate ion- 
ization balance in Ti I. Thus, inadequate results may reflect erro- 
neous photoionization cross-sections, which are computed in the 
hydrogenic approximation. We have performed a test decreasing 
the cross-sections for several low-excitation Ti I states, including 
a^¥, by three orders of magnitude. Although the individual line 
abundances decrease by ~ 0.03 dex, this is by far insufficient to 
bring two ionization stages in agreement. We conclude that it is a 
cumulative elfect of non-hydrogenic photoionization in Ti I, rather 
than cross-sections for individual atomic levels, that may lead to a 
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quantitatively different distribution of atomic level populations. In 
the absence of quantum-mechanical calculations, it is impossible to 
predict whether the hydrogenic approximation over- or underesti- 
mates the cross-sections for each Ti I level. Therefore, we can not 
explore this possibility further. 

If the problem is in NLTE, then it is likely that NLTE mod- 
elling of a similar atom, Fe, suffers from similar deficiencies as 
Ti. For example, accurate cross-sections for transitions caused by 
inelastic H I and e" collisions are missing for both atoms, and 
calibration of collision efficiency on observed spectra is neces- 
sary. In this case, one would expect the problem to be reduced 
or even eliminated if one compares the abundances determined 
from lines of equal ionization stages, i.e. Ti I and Fe I. The 
NLTE abundances based on Fe I lines and MAFAGS-OS model 
atmosphere^ are available only for HD 84937 and HD 122563, 
[Fei/H]= -2.00 + 0.07, respectively, [Fei/H]= -2.61 ± 0.09 (L. 
Mashonkina, private communication). The former value is 0.15 
dex higher than our standard Fe Il-based value for the metal-poor 
subdwarf, [Fe/H]= -2.15. For the metal-poor giant, the NLTE Fe 
abundance determined from Fe I lines is O.Il dex lower than the 
reference Fe Il-based value [Fe/H]= -2.50. It is interesting that 
Fe abundances derived from the NLTE analysis of Fe I lines are 
slightly discrepant with that derived from Fe II lines, and the sign 
of discrepancy is the same as for Ti both for the dwarf and giant 
stars. As a result, the [Ti/Fe] ratios for the two metal-poor stars de- 
termined from the NLTE analysis of Ti I and Fe I lines agree much 
better with that obtained from Ti II and Fe II lines. For a subdwarf 
HD 84937, we obtain [Tii/Fei]= -t-0.44 ± 0.02, which is consistent 
with [Tiii/Fen]= +0.35 ± 0.07 within the respective uncertainties 
of both values. For a giant HD 122563, [Ti,/Fe,]= +0.23 ± 0.15 
and [Tin/Fen] = +0.22 ± 0.02. Although the standard deviations of 
abundances are still large, the NLTE approach in this case has the 
advantage that it minimizes offsets between two ionization stages 
of Ti for both evolved and main sequence stars. 

There is also some similarity between our results, i.e. an ap- 
parent need for the high degree of th ermalization of Ti I levels i n 
the atmospheres of dwarfs, and that of lKom. Shi & GehrenI j2003h . 
who investigated kinetic equilibrium of Fe in metal-poor stars with 
the same methods, model atmospheres, and codes, as we do. They 
find that a very large scaling factor to Drawin's H I collision rates, 
5h = 3, in addition to the e nforced thermalization of the upper Fe 
I levels ( iGehren et alj200lh . is necessary to fuUfiU ionization equi- 
librium of Fe in their sample of dwarfs and subgiants with [Fe/H] 
< -I. Their results also indicate that NLTE effects on Fe II are 
negligible, which is supported by the results of Mashonkina et al. 
(2010) for Fe and is also true for Ti II (Table[6ll. 

We conclude that the ratios [Ti/Fe] in a stellar atmosphere 
computed exclusively using Ti II and Fe II lines should be robust, 
at least in the framework of ID modelling, and should be used 
for Galactic chemical evolution studies, whenever possible. It is 
strongly recommended to avoid Ti I lines, as well as to take av- 
erage of two ionization stages, in abundance analyses and in deter- 
mination of surface gravity for evolved stars. Otherwise, systematic 
offsets between giants and dwarfs are unavoidable and will produce 
spurious abundance trends with metallicity. 



iGehren et al] )2004l20o3) do not provide estimates of NLTE Fe abun- 
dances based on Fe I lines and MAFAGS-ODF model atmospheres for the 
stars in our sample. 



4 SUMMARY 

Statistical equilibrium of Ti was computed for a restricted range of 
stellar parameters, focussing on late-type stars. The Ti model atom 
was constructed with available experimental atomic data for levels 
and radiative transitions. Photoionization was assumed hydrogenic, 
except for the few Ti I levels, for wh ich we adopted experimental 
cross-sections from Yang et al. ( 2009) . Inelastic H I collision rates 
were compute d following the classicallDrawinI jl968l) formalism in 
the version of lSteenbock (&: Holwegej ll984h . The Drawin's cross- 
sections were increased by a factor of three, Sh = 3, because this 
value produced the smallest abundance scatter and satisfied ioniza- 
tion equilibrium of Ti I/Ti II for the Sun. 

NLTE effects on Ti I level populations and lines are signifi- 
cant for any combination of T^s, logg, and [Fe/H]. The dominant 
effect is overionization that makes NLTE particularly important at 
low metallicity. In the solar case, LTE approach underestimates the 
abundances by 0.05 - 0.1 dex. The solar Ti abundances derived 
from NLTE profile fitting of Ti I and Ti II lines with MAFAGS-OS 
theor etical model atmosphere are log e = 4.94 ± 0.05 dex (logg/'s 
from iBlackwell-Whitehead et ai]|2006h and log e = 4.95 ± 0.06 
dex (loggf's from IPickering et" I I2OOII) . Both values agree well 
with the Ti ab undanc e in C I meteorites, loge = 4.91 ± 0.03 dex 
jAsplund et alj|2009l) . The results obtained with ^/-values for Ti 
I transitions measured with Oxford furnace jBlackwell et al.|[l986l 
and earlier references) are identical, log e = 4.93 ± 0.04 d ex. The 
g/-values for Ti II transitions from ^ Bizzarri et al.l ( 1 19931) lead to 
~ 0.04 dex larger solar Ti abundances. MAFAGS-ODF models 
and/or LTE approach do not perform as good as MAFAGS-OS and 
NLTE with calibrated collision efficiency. 

Ti abundances were computed for four metal-poor stars: HD 
102200, HD 84937, HD 140283, HD 122563. LTE approach leads 
to an agreement between Ti I and Ti Il-based abundances for un- 
evolved stars, however it produces large discrepancy for giants. The 
same problem was also reported in other LTE studies of Ti abun- 
dances in metal-poor stars. Under NLTE, the discrepancy between 
the two ionization stages is significantly reduced for giants. How- 
ever, now the Ti I lines give higher abundances than the Ti II lines 
for dwarfs and subgiants. This failure should by no means be seen 
as invalidity of the NLTE approach, but rather as a manifestation 
of inaccurate or missing atomic data. It is probably important to in- 
clude highly-excited predicted levels and transitions between them 
in the model atom. In addition, accurate photoionization cross- 
sections and cross-sections for inelastic collisions with H I atoms 
are urgently needed. 

One of 'temporary' solutions might be to relate NLTE abun- 
dances determined from lines of equal ionization stages, i.e. Ti I 
and Fe I. In our preliminary tests, the [Ti/Fe] ratios for the metal- 
poor dwarf HD 84937 and giant HD 122563 determined from the 
NLTE analysis of Ti I and Fe I lines agree with those obtained from 
Ti II and Fe II lines. However, the standard deviations of abun- 
dances are not small. Also, one has to keep in mind the remaining 
discrepancy between Fe I and Fe II lines. 

At present, we strongly recommend to disregard Ti I lines in 
abundance analyses of giants, as well as in calibration of their sur- 
face gravities. Since the NLTE effects on Ti II and Fe II are small 
for the range of stellar parameters investigated in this work, the ra- 
tios [Ti/Fe] computed exclusively using Ti II and Fe II lines should 
be robust and should be used for Galactic chemical evolution stud- 
ies. 
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Table 6. Ti abundances for the selected sample. 



Object 


[Fe/H] 


[Mg/Fe] 




[Tii/Fe] ODF 




[Tiu/Fe 


;] ODF 


[Tii/Fe 


] OS 


[Tin/Fe] OS 










LTE NLTE 




LTE 


NLTE 


LTE 


NLTE 


LTE NLTE 


HD 84937 


-2.16 


0.32 


1 


0.49 ±0.02 0.63 ±0.02 


6 


0.34 ±0.08 


0.34 ±0.09 


0.46 ±0.02 


0.59 ±0.02 


0.35 ±0.06 0.35 ±0.07 


HD 102200 


-1.28 


0.34 


16 


0.28 ±0.03 0.39 ±0.05 


11 


0.28 ±0.05 


0.28 ±0.04 








HD 140283 


-2.38 


0.43 


5 


0.19 ±0.03 0.34 ±0.02 


6 


0.17 ±0.07 


0.18 ±0.08 








HD 122563 


-2.51 


0.45 


5 


-0.12 ±0.07 0.13 ±0.11 


3 


0.29 ±0.03 


0.34 ±0.03 


-0.19 ±0.6 


0.12±0.15 


0.21 ±0.02 0.22 ±0.02 
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